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ABSTRACT 



The high-mass X-ray binary (HMXB) SMC X-1 exhibits a superorbital mod- 
ulation with a dramatically varying period ranging between ~ 40d and ~ 60d. 
This research studies the time-frequency properties of the superorbital modula- 
tion of SMC X-1 based on the observations made by the All-Sky Monitor (ASM) 
onboard the Rossi X-ray Timing Explorer (RXTE).We analyzed the entire ASM 
database collected since 1996. The Hilbert-Huang Transform (HHT), developed 
for non- stationary and nonlinear time series analysis, was adopted to derive the 
instantaneous superorbital frequency. The resultant Hilbert spectrum is consis- 
tent with the dynamic power spectrum while it shows more detailed information 
in both the time and frequency domains. The RXTE observations manifest that 
the superorbital modulation period was mostly betweenn ~ 50 d and ~ 65 d, 
whenas it changed to ~ 45d around MJD 50,800 and MJD 54,000. Our anal- 
ysis further indicates that the instantaneous frequency changed in a time scale 
of hundreds of days between ~MJD 51,500 and ~MJD 53,500. Based on the 
instantaneous phase defined by HHT, we folded the ASM light curve to derive 
a superorbital profile, from which an asymmetric feature and a low state with 
barely any X-ray emissions (lasting for ~0.3 cycles) were observed. We also 
calculated the correlation between the mean period and the amplitude of the 
superorbital modulation. The result is similar to the recently discovered rela- 
tionship between the superorbital cycle length and the mean X-ray flux for Her 
X-1. 

Subject headings: accretion disks — stars: individual (SMC X-1) — X-rays: binaries 
—X-rays: individual (SMC X-1) 



-3- 



Introduction 



Many X-ray binaries exhibit long-term modulation with periods that are longer than 
their orbital periods. The modulation time scales of superorbital periods are typically 



tens to hundreds of days. A few X-ray binaries such as Her X-1 (JTananbaum. et al. 



I972I ) show superorbital modulation with stable periodicities, though re cent discoveries 



show that the superorbital cycle length may vary within a small range ( jLeahy fc Igna 
20101). A small group of X-ray binaries show quasi-periodic superorbital modulation 
with periods changing w ithin a certain range. The most prominent of these is SMC X-1 



( IClarkson et al. 



2003a| ). Some other X-ray bi naries such as Cyg X-2 exhibit superorbital 



modulation with multi-periodic modulations fjClarkson et al 



2003bl ). Such long-term 



superorbital variations rnay be caused by the pre cession of a tilted or warped accretion disk 



( IWijers fc Pringle 



1820-30 ( 



1999 



Ogilvie fc Dubus 



Chou &: Grindlay 



2OOII ). a hierarchical third companion like 4U 



200ll ). or by reasons currently unknown. 



The high-mass X-ray binary (HMXB) SMC X-1, w hich was first discovered in 1971 



f lLeong et al. 



19711 ) . consists of a 1.06 Mq ne utron star 



BO I type supergiant with a mass of 17.2 Mq ( iReyno 



van der Meer et al. 



the neutron star is 0.71 s and it decreases with time (IWoidowski et al. 



ds et al. 



1993) 



20071 ) and a 



period of this system is ~ 3.89 d as evaluated by its ec 



orbital incli nation angle is ~ 67° ( 



time delay. 

From t 
(HEAO-1), 



Wojdowski et al. 



van der Meer et al. 



ipse 



he sp in period of 



( ISchreier et al. 



19981). T he orbital 



I972I), and the 



20071 ). By analyzing the mid-eclipse 



(119981 ) found that the orbital period also changes with time 



-3.35 X lO'V"^) 



le observational data c ollect ed by the High Energy Astronomy Observatory 1 



Gruber fc Rothschild 



a rough time scale of 60 d. 



(119711) found that the X-ray flux of SMC X-1 varies with 



Wojdowski et al. 



( 119981 ) confirmed this periodicity and showed 



that the superorbital cycle length was not stable during the first ~ 500 d observations made 
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by the All-Sky Monitor (ASM) onboard the Rossi X-ray Timing Explorer (RXTE). After 



the ASM collected sufficie nt amounts of dat a, several time-frequency ana 



as the wavelet transform (JRibo et al 



2OOII ). d ynamic power spectrum 



ysis methods su ch 



(jClarkson et al. 



2003al ). and sliding Lomb-Scargle periodogram (JTrowbridge et al. 1120071 ) were applied to 
the light curve to investigate the variations in the superorbital period of SMC X-1. The 
mechanism of superorbital modulations in SMC X-1, similar to those in Her X-1, LMC 
X-4, and Cyg X-2, are interpreted by a warped and tilted accretion disk ( IWojdowski et al. 
19981 ). When the disk precesses, it obscures our line of view to the central X-ray source. 
However, the mechanisms that cause the precession period to change with time are still 
unknown. 



Ribo et al. 



(120011 ) analyzed the first ~ 1700 d of the ASM light curve and found that 
the superorbital period first decreases from ~ 60 d to ~ 45 d and subsequently increases to 
~ 60 d. By fitting the maximum points of the wavelet spectrum with a si nusoidal curve. 



they fo und that the modulation period varies with a period of 1421 ± 8 d. 



Clarkson et al. 



(|2003aJ) analyzed the first ~ 2200 d of the ASM light curve along with all the Burst and 
Transient Source Experiment (BATSE) data using the dynamic power spectral technique. 
The dynamic power spectra of both the RXTE and BATSE light curves show that the 
superorbital modulation period changes between ~ 60 d and ~ 45 d. The superorbital 



Trowbridge et al. 



(120071 ) analyzed 



period appears to vary in a ~ 1600 d time scale, 
the first ~4000 d of ASM data using the sliding Lomb-Scargle periodogram and the the 
cycle length counting method, and the results indicated that the cycle length variation is 
perhaps not as smooth as that obtained from the dynamic power spectrum. The cycle 
length variation possibly contains a shorter periodicity in addition to the ~ 1600 d period. 
However, due to the constraints of the window size in both the dynamic power spectrum 
and the sliding Lomb-Scargle periodogram, it is difficult to further analyze the periodicities 
of the superorbital period variation. 
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In order to study the superorbital period variatio n with sufficient sta tistical accuracy, 



19981 ) to derive 



here we adopt the Hilbert-Huang transform (HHT) (jHuang et al. 
instantaneous frequency of the ASM light curve data for periodicity analysis. HHT is a 
new time-frequency analysis method, designed for non-stationary and nonlinear time series 
analysis. It consists of the empirical mode decomposition (EMD) and the Hilbert spectral 
analysis. In this work, we first apply EMD to decompose the ASM light curve into a number 
of intrinsic mode functions (IMFs) which are suitable for the Hilbert transform, and then 
calculate the instantaneous frequency and amplitude of certain IMFs, which form the basis 
of further analysis on periodicity of superorbital period variation, folded light curve, and 
correlations between the period and the amplitude. The proposed analysis allows us to 
explore subtle variations in the superorbital periods which were not detected by existing 
methods. 

The rest of this paper is organized as follows. In section |2l we briefly introduce the 
ASM observations and the algorithm of HHT. Section [3] describes the time-frequency 
analysis of ASM observations based on HHT. Finally, we discuss our results in section HI 

2. Observation and Data Analysis 

2.1. RXTE ASM 

Since the RXTE was launched in late 1995, the ASM continuously sweeps the entire 
sky once every 90 minutes. The energy range of the ASM is 1.3 to 12.1 keV. This can be 
further divided into three energy channels (chl: 1.3-3.0 keV, ch2: 3.0-5.0 keV, and ch3: 
5.0-12.1 keV). The summed band data collected from 1996 to 2010 (MJD 50134 to MJD 
55371), with a total time span of ~ 5000 d were used in this time-frequency analysis. The 
dwell data, where all the eclipses were removed according to the ephemeris proposed by 
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Wojdowski et al. 



method of 



1998 



), we r e binn ed into a one-day averaged light curve. Following the 



Trowbridge et al. 



( 120071 ). we did not exclude those data points that have a count 
rate falling below zero. Figure [T] shows the ASM light curve of SMC X-1. Although the 
dwell data were binned into a one-day averaged light curve, it still contained some missing 
points or small observ ation gaps. We then used the piecewise cubic Hermite interpolation 



( jKahaner et al. 



19891 ) to interpolate the ASM light curve to evenly sampled data. 



2.2. Hilbert-Huang Transform 

Owing to the high variability of the superorbital cycle length of SMC X-1, it is suitable 
to analyze the time-frequency properties of its superorbital modulation by HHT. We use 
the Hilbert transform to calculate instantaneous frequency of the data. For a time series 
x{t), this is achieved by first calculating its conjugate yit), defined as 



yit) 



-p 



11 






dt'. 



where P indicates the Cauchy principal value. In the complex plane, we define an analytic 
signal, z{t), as 

z{t) = x{t) + iy{t) = a{t)e''^'\ (2) 

where a{t) = ^/x^tf+y^t) is the amplitude and 9{t) is the angular phase function. The 
instantaneous angular frequency u{t) is then defined as the time derivative of 0{t), 

de{t) 



uj{t) 



dt 



(3) 



Accordingly, we can obtain the frequency instantaneously as long as the sampling rate is 
sufficiently high . However, the insta ntaneous frequency is meaningful only if the time series 



x{t) is an IMF (JHuang et al. Ill998l ). An IMF satisfies the following two conditions: (1) the 



number of extrema and the number of zero crossings must be identical or differ by one; and 
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(2) at any of the points in the series, the mean value of the upper envelope ( defined by the 



loca 



maxima) and the lower envelope (defined by the local minima) is zero f lHuang et al 



19981 ). Because the real data are usually not IMF, it is difficult to directly apply t he Hi lbert 



Huang et al. 



fll998h 



transform on primary data in most cases. To address this issue, 

proposed the EMD method to decompose a time series into several IMFs such that the 

Hilbert Transform is applicable. 

The EMD method assumes that any time series consists of some oscillatory 
components. The decomposition scheme utilizes the actual time series for the construction 
of the decomposition base rather than decomposing it into a prescribed set of base functions. 
The decomposition is achieved by iterative "sifting" processes for extracting mod es by 



identification of local extrema and subtraction of local means ( jHuang et al. 



1998 



). The 



iterations are te rminated by a criteri on of convergence. For details of the sifting, reference 



is here made to 



Huang et al. 



(119981 ). Under the procedures of EMD, the original time 



series x{t) is decomposed according to 



xit) = y^Q + rn, 



(4) 



i=l 



where q's are IMFs and r„, is a residue. For a time series with infinite length, the 
decomposition satisfies the properties that the decomposed components are orthogonal to 
each other and form a complete set. For a time series with finite length; however, the 
orthogonality of the decomposition may not strictly hold. In this case, we minimize the 
orthogonality index (01) defined by 

EtEi^jCi(t)cj-(t) 



01 



Et^'it) 



(5) 



to make the decomposition "unique" . Furthermore, if the original data is noisy or contains 
high-frequency intermittent signals, the decomposition may suffer the mode mixing problem , 



i.e. the modulation signal with the same time scale across different IMFs (JHuang et al. 



19981 ). This would cause many spurious, confused signals in the transition region. When 
the mode mixing problem occurs, it is difficult to obtain a rigorous instantaneous frequency, 
ui(t), of the modulatio n to be further analyzed. Then, the ensemble EMD (EEMD) 



(IWu and Huang 



20091 ). a revised version of EMD implemented by taking an ensemble 



average of the components decomposed from combinations of primary time series and 



distinc t white noises, is used to mitigate the mode mixing problem, and a post-processing 

m 



EMD flWu and Huang 



20091 ) is applied afterward on the outcome of EEMD to guarantee 



that the final IMFs fulffi orthogonality and completeness. 

After decomposing the time series into IMFs, the Hilbert transform is applied on each 
IMF. As a result, the data can be expressed as 



x{t) = y^ CLjit) exp I i I 

.7 = 1 ^ "^0 



ujj{t')dt' 



(6) 



where aj{t) and (jJjit) are the instantaneous amphtudes and frequencies of decomposed 
components. The amplitude and frequency of each component in Eq.([6]) are functions of 
time, instead of constant amplitude and constant frequency in the Fourier transform. 

Figure |5] shows the decomposition of the interpolated ASM light curve of SMC 
X-1. There are 11 IM Fs, 1 residue component, and 01 = 0.036. The significance test 



( IWu and Huang 



2004 ) on these IMFs suggests that the 40-60 d superorbital modulation 



signal of SMC X-1 mostly concentrates in the fifth component (cs). The values of all the 
other components are lower than or marginal at the 3cr white noise level, suggesting that 
they act in a manner similar to white noise. These components are probably caused by the 
observational noise associated with non-periodic modulations in other time scales. 
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Results 



3.1. Hilbert Spectrum 



After decomposin g the ASM light curve into IMFs, the normahzed Hilbert transform 



(JHuang fc Long 



20031 ) was applied on the IMFs to obtain the instantaneous frequencies 



and amplitudes. The normalized Hilbert transform was proposed to overcome th e limitation 



of ab sence of strong amplitude modulations imposed by the Bedrosian theorem (JBedrosian 
19631 ). The resultant Hilbert spectrum is a three-dimensional map which displays how the 
modulation period and amplitude vary with time. Figure E] shows the result, in which the 
frequency range is divided into 3000 bins and the spectrum is smoothed by a Gaussian 
filter for clarity. The color map represents the Hilbert energy spectrum with the magnitude 
of energy defined as squar e of the amplitude. For comparison, the dynamic power spectral 



technique, as described in 



Clarkson et al. 



(J2003al ). was also applied on the ASM light 



19821 ) of the first 



curve. We first calculated the Lomb-Scargle power spectrum (iScargle 

200 d of the ASM data, then successively moved the 200 d data window forward by a step 

size of 10 d and applied the Lomb-Scargle algorithm on each segment to obtain a series of 



power spectra. The power spectra are plotted in Figure [3] i n the contour 



spectral power as defined in the Lomb-Scargle periodogram (IScargle 



brm with the 



19821 ). It is obvious 



that the Hilbert energy spectrum is consistent with the dynamic power spectrum. Because 
the spectral analysis based on HHT is independent of window size, more detailed structures 
in the light curve data can be observed from the high-resolution Hilbert spectrum than the 
dynamic power spectrum. 

The first remarkable feature of both spectra is that the superorbital modulation period 
changes dramatically from ~ 60 d to ~ 45 d and goes back to ~ 60 d bet ween ~MJD 



50,200 and ^MJD 51 , 100. This feature is consistent with previous studies (JRibo et al. 



2001 



Clarkson et al. 



2003a; 



Trowbridge et al. 



20071 ). Interestingly, a similar variation in 



-10- 

the superorbital modulation period is repeated between ~MJD 53,600 and ~MJD 54,400. 
The time durations of these two events are roughly 800 d and the separation between them 
is ~ 3200 d. This indicates that such a phenomenon is probably recurrent. In addition, 
although the superorbital period is relatively stable, it seems to oscillate between P ~ 50 d 
and ~ 65 d over a time scale of hundreds of days between ~MJD 52,000 and ~MJD 53,000. 
This feature can be seen very marginally in the dynamic power spectrum, but it is difficult 
to identify due to the limitation of the window size. By contrast, the Hilbert spectrum 
based on HHT displays the instantaneous frequency (or period) more accurately, enabling 
the application of further timing analysis methods on the instantaneous frequency to assess 
the periodicity hidden in the superorbital period change. 

In order to study the periodicities in the variation of the superorbital frequency, the 
Lomb-Scargle periodogram was applied on the instantaneous frequency of the component 
C5, and the result is shown in Figure HI The highest p eak is located at P = 1674 d; this 



is con sistent with the previously reported periodicity (JRibo et al. 



2001 



Clarkson et al. 



2003al ). However, this peak is the second harmonic of P ~ 3200 d, and it corresponds to 
the separation of the events occurring when the superorbital period changes between ~ 60 
d and ~ 45 d. Moreover, a few peaks with relatively less power values located at P ~ 240 d 
and P ~ 320 d are observed. These peaks are likely due to the oscillation feature mentioned 
in the previous paragraph. 

It is possible to further apply time-frequency analysis methods such as HHT and the 
dynamic power spectrum on the instantaneous frequency to examine the variation of the 
superorbital modulation with time. Figure [5] shows a dynamic power spectrum with a 
window size of 1,000 d and a movement of 10 d. The most prominent peaks are located at 
around MJD 50,800 and MJD 54,000. These peaks correspond to the short-period state 
of the superorbital modulation in the time domain in Figure |3] and the highest peaks in 
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the frequency domain in Figure HI During ~ MJD 52,000 and ~ MJD 53,500, the most 
prominent periodicities disappear and smaller periodicities from P ~ 240 d to P ~ 320 d 
are observed. These correspond to the oscillation feature of the instantaneous frequency 
with a time scale of hundreds of days, as shown in Figure [31 Figure [5] shows that the 
oscillation feature seems to appear again after ~ MJD 54,500, indicating that it might exist 
in the entire observation time span, except during the two events where the period changes 
dramatically between ~ 60 d and ~ 45 d. 

3.2. Superorbital Modulation Profile 

It is difficult to derive a proper profile to describe the characteristics of superorbital 
modulation for highly v ariable cycle lengths from a light curve folded with a fixed period. 



Clarkson et al. 



( l2003al ) first selected a segment of light curve at which the frequency seems 
stable in the dynamic power spectrum and subsequently folded it with a fixed period. 
However, the period still changes during the entire observation time even in the stable region 



of the dynamic power s pectrum, according to the Hilbert spectrum anc 



cycle length counting by 



Trowbridge et al. 



( 120071 ). As an improvement. 



the analysis of 



Trowbridge et al. 



(120071 ) fitted the smoothed light curve with fourth-order polynomial functions to obtain the 
minimum points of the superorbital modulation, and folded the whole light curve according 
to the ephemeris defined by minima. The difficulty then becomes the identification of some 
minima for the presence of noise in the ASM light curve and distortions of the modulation 
shape when the period changes dramatically. By contrast, the light curve folded according 
to the phase defined as 0(t) = 6(t)/27T in HHT does not suffer the above-mentioned 
problems. Figure El shows the folded light curve of SMC X-1, and the phase zero epoch is 
defined by the first data point in the light curve (MJD 50,134). 

The folded light curve shown in Figure [61 seems mostly in agreement with the results 



- 12- 

of iTrowbridge et al. I (l2007l ). The major distinction between this superorbital modulation 
profile and that obtained from previous studies is a clear low state with a negligible X-ray 
flux lasting for ~ 0.3 cycle. Moreover, the asymmetric features of the superorbital profile 
can also be obtained in the folded light curve. In order to describe the asymmetric profile, 
the low state count rate was defined as the mean count rate during phase ~ 0.2 to ~ 0.55 
and the high state count rate was considered by averaging the count rate between phase 
~ 0.75 to ~ 1.0. Thus, the rising and falling time scale can be estimated by calculating the 
time scale between 10% and 90% of amplitude. By this definition, the rising time scale is 
~ 0.14 cycle (from phase 0.05 to 0.19), which is shorter than the falling time scale of ~ 0.19 
cycle (from phase 0.54 to 0.73). 

The ASM light curve can further be divided into three channels. The spectral hardness 
provides crude information about the emission mechanisms during superorbital modulation. 
The ASM hardness ratio was defined as ch3/{chl + ch2), since ch3 and chl + ch2 have 
similar photon count rates during the high state. Only the data with a signal to noise ratio 
(SNR) greater than 5 were chosen and folded according to the superorbital phase defined in 
HHT. During the entire observation time span, there are no data points with SNRs greater 
than 5 between phase 0.73 and 1.05. Therefore, no reliable hardness ratio in the low state 
was obtained. The folded hardness ratio is shown in Figure [71 It is easily observed that 
the hardness ratios in the high state (from phase 0.19 to 0.54) are relatively stable. In 
the transition state (phase 0.05-0.19, and 0.54-0.73), the hardness ratios do not appear to 
significantly deviate from the mean value as compared with that in the high state although 
the errors are much greater. 
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3.3. Correlation between Period and Amplitude 

Another remarkable feature observed in this study is the correlation of the amplitude 
with the period. Using the instantaneous frequency and amplitude derived by HHT, the 
correlation between the cycle-averaged period and the Hilbert amplitude in the most 
significant IMF (cs) was calculated. The rank correlation coefficient is r = 0.46 with a null 
hypothesis probability value of 3.2 x 10~^, indicating a significant correlation between the 
period and the amplitude (see Figure [H]). 

In order to obtain the correlation between the period and amplitude in the original 
data set, the rms amplitude of each cycle was calucalted. The result shows that the period 
and the rms amplitude still show a significant correlation (r = 0.32 with null hypothesis 
probability value of 1.4 x 10"'^, see Figure [9]). Although the correlation appears to change 
to an anti-correlation after P > 60 d, the number of cycles is too few to obtain a significant 
correlation for periods l onger than ~ 60 d. Similar correlation is observed in the superorbital 



modulation of Her X-1 (JLeahy fc Igna 



2010 



Still fc Bovd 



2004h . 



4. Discussion 

We have successfully performed HHT-based time-frequency analysis on the ASM light 
curve of SMC X-1. The Hilbert spectrum manifests variations in the instantaneous period 
of the superorbital modulation. The instantaneous phase defined by HHT allows us to fold 
the light curve to derive a reasonable superorbital profile. The time-frequency information 
obtained from the HHT analysis is more abundant than the traditional Fourier-based 
methods. This benefits our investigations on the superorbital modulation of SMC X-1. 

Recent observations by the Chandra and XMM- Newton show that the emission of SMC 
X-1 could be described as a blackbody component with kTsB ~ 0.18 keV and a hard power 



-14- 



law component (JNeilsen et al. II2004I ). The emissions with energy < 1.0 keV are dominated 
by the blackbody radiation. These emissions are reprocessed X-rays emerging from the 



illuminated warp region of the accretion disk fJHickox fc Vrtilek 



20051 ). On the other hand, 



the high energy emissions are dominated by the power law components owing to pulsar 
beams emerging from the neutron star surface. Therefore, the energy range of all the ASM 
channels is dominated by the power law components. From the folded hardness ratio, we 
could not obtain the absorption feature since the hardness ratios in the transition state do 
not significantly differ from that in the high state. This indicates that the warp region of 
the accretion disk is optically thick for all the ASM energy bands, and the superorbital 
profile is determined by the covering of the central source by the warp region. 

Furthermore, we found a significant correlation between the amplitude and the period 
of the superorbital modulation. A similar correlation is observed in Her X-1, which is a 
Low-Mass X-ray Binary (LMXB) with an o rbital period of 1.7 d and s hows a relatively stable 



superorbital modulation in a 3 5 d period (ITananbaum. et al 



obtained by the O— C method. 



19721 ). From the ephemeris 



Still fc Boydl ( 12004J ) showed that the time-averaged peak 



main-on flux and the precess ion period have a pos itive correlation though the quantity of 



obtained data is very small. 



Leahy &: Igna 



( I2OIOI ) analyzed the ASM light curve using 



the cross-correlation method and showed that the cycle length and average ASM count 
rate indeed have a significant correlation. Althou gh the mechanisra s that cause the disk 



Wijers fc Pringle 



f ll999h showed that for 



precession period change with time are unknown 
a radiation-induced warp disk, the precession period of the outside edge of the disk would 
anti-correlate with the strength of the radiation field. From the observation that the low 
state of SMC X-1 has barely X-ray flux, the modulation amplitude might be associated 
with the central X-ray flux if the central emission source is totally unc overed during the 



high state. Our resu lts for SMC X-1, as well as the results for Her X-1 flStill fc Bovd 



Leahy fc Igna 



2OIOI ). seem to be inconsistent with this prediction. 



Leahy fc Igna 



2004; 



teoiol 
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accounted for the average ASM count rate, superorbital cycle length, and spin-down of 
the neutron star in their studies and obtained a consistent result, indicating that the cycle 
length could directly be related with the variations in accretion flow. Fo r SMC X-1, the 
evolution of the pulse frequency was demonstrated by 



Inam et al. 



(120 lOl ). The frequency 



residual during the RXTE observations seems to vary with the time scale of the superorbital 
variation. However, the observed sampling is insufficient to obtain a significant correlation 
between the spin frequency residual and the superorbital period variation. 

An alternative possibility is that the amplitude might be affected not only by the 
X-ray flux of the c entral source but also by the variations in tilt angle, warp mode, etc. 



Ogilvie &: Dubus 



(120011 ) studied the stability of radiation-driven warping of accretion 
disks in X-ray binaries and found that the binary separation (r^) and mass ratio (g) are 
good indicators of warping behavior. SMC X-1 lies in the marginal region of the r^ — q 



paramet 



er space where mode stable warping and other warping modes can possibly 



appear (jClarkson et al. 



2003bh . In contrast to SMC X-1, Her X-1 and LMC X-4 exhibit 



only stable mode warping. However, the stable warping in Her X-1 also shows cycle length 
variations. This indicates that mode is perhaps not stable, and the correlation between 
cycle length and amplitude or mean count rate is a property of mode 0. 



This research has made use of data provided by the ASM/RXTE teams at MIT and 
at the RXTE SOP and GOF at NASA's GSFC. This research is supported by grant NSC 
97-2112-M-008-012-MY3 of the National Science Council of Taiwan. 
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Fig. 1.— All-Sky Monitor (ASM) light curve of SMC X-1. 
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Fig. 2. — Intrinsic mode functions (IMFs) decomposed by ensemble empirical mode decom- 
position (EEMD). The fifth IMF (cs) is responsible for the ~ 40 d to ~ 60 d superorbital 
modulation. 
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Fig. 3. — Hilbert energy spectrum (color map) and dynamic power spectrum (contour plot). 
The blue curve represents the instantaneous frequency of IMF c^ after Gaussian smoothing 
and the color depth denotes the magnitude of the Hilbert energy. The window size of dynamic 
power spectrum is 200 d and the movement is 10 d. 
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Fig. 4. — Lomb-Scargle periodogram of the instantaneous frequency of IMF c^. The highest 
peak is located at / ~ 0.0006 (or P = 1647 d), which is the second harmonic of P ~ 3200 
d peak. The other prominent peaks are located at / ~ 0.0031 (P ~ 240 d) and / ~ 0.0031 
(P ~ 320 d). 
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Fig. 5. — Dynamic power spectrum of the instantaneous frequency of IMF C5. 
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Fig. 6. — Folded light curve of SMC X-1 according to the phase defined by HHT. The number 
of bins in one cycle is 40. 
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Fig. 7. — Hardness ratios folded according to the superorbital phase. Only the data points 
with a signal to noise ratio greater than 5 were chosen to calculate the hardness ratios. 
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Fig. 8. — Correlation of cycle-averaged period and Hilbert amplitude of IMF C5. 
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Fig. 9. — Correlation of cycle-averaged period and rms amplitude of original ASM light 
curve. 



